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I.  INTRODUCTION 


Figure  1.  Controlling  the  Behavior  of  a  Chaotic  System 


The  objective  of  a  ballistic  missile  defense  (BMD)  sensor  system  is  to  find  enemy  targets, 
aim  weapons,  and  determine  whether  the  weapon  destroyed  the  selected  target.  This 
objective  is  broken  into  three  respective  functions:  discrimination,  pointing  and  tracking, 
and  kill  assessment.  Sensitive  optical  and  electronic  equipment  aboard  interceptors  and 
satellites  demand  aiming  accuracies  on  the  order  of  a  tenth  of  a  micro  radian.  Vibrations 
caused  by  thermo-mechanical  flutter,  thruster  firing,  and/or  structure  borne  noise 
significantly  reduce  sensor  precision.  Suppression  systems  are  therefore  vital  for  the 
development  of  highly  accurate  BMD  weapon  systems. 

Our  intention  is  to  develop  an  innovative  control  vibration  device  called  a  Chaotic  Motion 
Stabilizer  for  BMD  applications.  Chaotic  Motion  Stabilizers  transform  chaotic  behavior 
into  predictable  periodic  motion  as  Figure  1  illustrates.  Our  product  enhances  the 
performance  of  active  damping  systems,  eliminates  the  need  for  damping  devices  in  some 
applications,  and  increases  the  variability  of  passive  damping  methods. 

Chaotic  Motion  Stabilizers  could  make  active  damping  systems  a  more  attractive 
alternative  in  some  applications,  and  improve  the  effectiveness  and  simplify  the  control 
algorithm  in  others.  Active  damping  uses  anti-noise  to  cancel  extraneous  vibrations.  This 
entails  adding  anti-noise  or  noise  180  degrees  out  of  phase  with  the  unwanted  vibrations 
to  mirror  the  incoming  noise.  When  the  two  add  together,  the  result  is  silence.  There  are 
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two  active  approaches  —  feedfoward  and  feedback  control  systems.  In  both  of  these 
systems,  the  effectiveness  of  the  technique  hinges  on  the  systems’  ability  to  predict  the 
incoming  noise  at  the  location  where  anti-noise  is  applied. 

In  a  feedback  scheme,  a  down  stream  microphone  returns  a  sampled  signal  to  the 
controller.  The  controller  then  sends  an  anti-noise  signal  to  the  speaker  that  neutralizes 
the  incoming  vibration.  A  feedback  system  will  only  work  for  predictable  periodic 
motions.  The  fundamental  feedback  problem  is  the  inability  to  predict  the  chaotic 
behavior  of  an  upstream  location.  Figure  2  clearly  illustrates  this  problem.  A  paper 
stream  with  randomly  colored  black  and  white  dots  is  moving  at  a  constant  rate  across  a 
plane.  If  we  know  the  color  of  the  dot  at  location  1,  it  tells  us  nothing  about  the  color  of 
the  dot  at  location  0,  upstream,  because  the  system  is  random.  The  only  way  to  know  the 
behavior  of  an  upstream  location  is  when  the  dots  are  periodically  placed.  In  the  same 
manner,  active  feedback  can  only  work  in  situations  where  noise  is  periodic  or  predictable. 
When  faced  with  chaotic  motions,  a  feedfoward  system  is  more  appropriate. 

In  feedfoward  active  control  an  upstream  microphone  samples  the  magnitude  and  phase  of 
the  incoming  noise.  A  controller  uses  the  sample  to  drive  a  downstream  speaker  that 
provides  the  anti-noise.  A  feedforward  systems  effectiveness  is  dependent  on  the  ability  to 
predict  the  magnitude  and  phase  of  the  noise  at  the  downstream  speaker  location.  The 
method  works  for  both  chaotic  and  periodic  motions  assuming  two  favorable 
circumstances.  First,  they  can  create  a  time  delay  from  the  sampled  noise  to  the  speaker 
and  the  time  delay  is  slow  enough  for  the  control  electronics  to  act.  Second,  the  sampled 
motion  upstream  does  not  change  before  it  reaches  the  anti-noise  speaker.  Figure  3  more 
clearly  illustrates  feedforward  control.  The  same  paper  stream  in  Figure  2  is  moving  at  a 
constant  rate  of  1  cm/sec  across  a  plane.  If  we  know  the  dot  is  black  at  location  1,  in  1 
second  the  dot  will  still  be  black  at  location  2,  1cm  downstream.  We  can  use  the  same 
analysis  for  feedforward  active  control.  Feedforward  systems  have  shown  success  in 
neutralizing  chaotic  noise  in  air  ducts  and  mufflers. 

Feedback  and  feedforward  active  systems  are  not  effective  in  many  applications.  Active 
feedback  systems  are  not  capable  of  neutralizing  chaotic  noises  that  appear  in  most 
vibration  applications.  They  can  only  cancel  periodic  signals.  Feedforward  systems  can 
neutralize  chaotic  motions,  but  the  necessary  delay  times  required  for  effective 
performance  are  too  short  in  structural  applications.  The  wave  speeds  in  solids  are  faster 
than  the  electronic  algorithms  needed  for  digital  analysis  -  the  unwanted  vibrations  pass 
the  antinoise  actuator  before  it  has  a  chance  to  analyze  the  upstream  sensor  sample.  In 
addition,  many  structural  applications  cannot  use  one-dimensional  noise  source 
approximations  as  feedforward  systems  depend  upon.  At  a  single  location  on  a  structure 
the  source  of  noise  is  three  dimensional  and  could  come  from  a  combination  of  the  infinite 
number  of  directional  vectors  surrounding  that  location.  The  inherent  nature  of  chaotic 
control  is  not  dependent  on  a  noise  source,  in  fact,  it  relies  on  the  complexities  of  the 
incoming  noise.  With  their  ability  to  transform  chaotic  behavior  to  predictable  periodic 
motion,  Chaotic  Motion  Stabilizers  make  active  feedback  possible  eliminating  the  need  for 
feedforward  systems  altogether. 
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Figure  2.  Feedback  Phenomenon  of  a)Random  Dots  b)Periodic  Dots 
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Chaotic  Motion  Stabilizers  could  also  replace  damping  systems  in  some  electronics  and 
instrumentation.  Employing  a  stabilizer  could  eliminate  the  need  for  such  devices  since  on 
board  electronic  instrumentation  could  compensate  for  the  predictable  motions.  End  users 
would  prefer  semi-active  electronic  compensation  versus  active  cancellation  because  of  an 
active  system’s  adverse  collateral  effects.  These  include  introducing  vibrations  into  other 
areas  of  the  system  and  increasing  the  vibration  to  catastrophic  levels  if  the  system  fails.  If 
a  Chaotic  Motion  Stabilizer  fails,  it  is  acts  as  a  passive  damper. 

They  also  offer  an  innovative  approach  to  passive  damping  approaches  because  of  their 
variability.  The  passive  damping  approach  requires  the  application  of  vibration  absorbing 
materials  either  housed  in  a  damper  or  directly  applied  to,  the  desired  system.  Vibration 
absorbent  materials  include  fiberglass,  acoustic  foams,  rubber,  polymer  blends,  and  fluids. 
When  strained  these  materials  absorb  the  mechanical  energy  and  transform  it  into  heat; 
thereby  reducing  some  of  the  vibrational  response.  The  limitation  of  the  method  rests 
with  the  damping  material’s  mechanical  behavior.  Their  damping  properties  are 
dependent  on  temperature  and  frequency  where  optimal  damping  occurs  at  single 
location.  With  a  Chaotic  Motion  Stabilizer,  one  can  select  infinite  numbers  of  periodic 
orbits  that  could  improve  system  performance  in  a  diverse  range  of  environments. 

Table  1  summarizes  the  strengths  and  weaknesses  of  the  different  approaches  to 
vibrational  damping. 


Enhances  the  performance  of  active 

Unproven  technology 

systems 

•* 

Costly 

Chaotic  Motion  Stabilizers 

Variable  damping  through  selection 
of  infinite  numbers  of  periodic 
orbits 

Uses  a  semi-active  approach  that 
won’t  cause  catastrophic  failure 

Low  Cost 

Limited  range  of  operational 
effectiveness 

Passive  Damping  Devices 

Simple  to  apply 

Ineffective  for  many  applications 

Eliminates  unwanted  predictable 

Cannot  neutralize  chaotic  motions 

Active  Damping  Systems 

vibrations  completely 

in  many  applications 

Costly 

Table  1.  Competitive  Matrix  of  Damping  Systems 
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IL  METHODOLOGY 


Figure  4.  Dynamic  Behavior  of  an  Undamped  Swinging  Pendulum 


In  the  past,  scientists  believed  they  could  predict  the  future  indefinitely  by  knowing  the 
laws  of  nature  and  the  initial  conditions  of  the  system.  For  example,  with  a  knowledge  of 
the  laws  of  gravity,  Newton’s  laws,  and  the  planets  initial  conditions  —  position,  mass,  and 
velocity;  astronomers  can  predict  planetary  motions  thousand  of  years  into  the  future.  The 
inability  to  predict  the  future  in  other  systems  such  as  the  weather  and  turbulence  was 
blamed  on  a  lack  of  knowledge  of  the  laws  governing  the  system,  the  complexity  of  the 
system,  and/or  the  inability  to  measure  initial  conditions.  They  assumed  that  if  we 
obtained  all  of  this  information,  one  could  predict  the  behavior  of  seemingly  random 
systems  much  like  planetary  motion. 

Chaos  theory  killed  this  idea.  Simple  deterministic  systems  with  only  a  few  elements 
generate  seemingly  random  behavior.  We  now  accept  that  unpredictability  as  fundamental 
to  chaotic  systems.  The  reason  --  a  chaotic  systems’  sensitivity  to  initial  conditions.  Two 
identical  systems  with  very  small  differences  in  initial  conditions  diverge  exponentially. 
Minute  differences  become  so  amplified  that  although  the  system  is  predictable  in  the  short 
term,  its  long-term  behavior  is  completely  unpredictable. 

State-space  representations  are  useful  tools  for  observing  and  explaining  chaotic  systems. 
It  is  an  abstract  space  whose  coordinates  are  governed  by  the  number  of  initial  conditions 
needed  to  solve  the  differential  equations  modeling  the  system.  The  swinging  pendulum 
illustrates  a  simple  example  of  state-space.  It  has  two  state  variables:  angular  position 
and  velocity.  In  the  undamped  case,  Figure  4,  the  pendulum  swings  back  and  forth  for  all 
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Figure  5.  State-space  Representation  of  a  Swinging  Pendulum 


of  time.  The  position  and  velocity  states  with  respect  to  time  are  shown.  A  state-space 
representation  combines  the  two  as  Figure  5  illustrates.  The  steady-state  response 
represents  a  linear  attractor  whose  trajectories  are  bound  within  a  region  of  space,  in  this 
case  a  circle.  The  state-space  representation  of  a  chaotic  attractor  is  shown  in  Figure  6. 
The  axes  represent  states  of  the  system  and  the  loops  represent  the  dynamic  trajectory  of 
the  states  with  respect  to  time.  The  paths  never  retrace,  but  are  bound  within  regions  of 
state-space  known  as  chaotic  attractors.  These  attractors  depict  collections  of  many 
orderly  behaviors,  none  of  which  exclusively  dominates. 

We  have  based  our  product  on  the  observance  that  chaotic  attractors  typically  Jiave  infinite 
numbers  of  unstable  periodic  orbits.  The  device  attempts  to  stabilize  existing  orbits  in  the 
attractor  —  not  create  new  ones.  The  process  is  straight  forward.  First,  determine 
unstable  periodic  orbits  embedded  in  the  chaotic  attractor.  Examine  these  orbits  and 
choose  ones  that  yield  improved  system  performance.  Finally,  tailor  small-time  dependent 
parameter  perturbations  to  stabilize  the  orbit. 

To  make  the  necessary  system  parameter  changes  that  transform  chaotic  motion  into 
periodic  motion,  we  will  utilize  electrorheological  (ER)  materials.  These  materials  are 
suspensions  of  dielectric  particles  in  non-polar  liquids  that  exhibit  dramatic  reversible 
changes  in  rheological  response  when  exposed  to  an  electric  field.  A  simplified  physical 
description  of  the  behavior  is  that  the  material  transforms  from  a  liquid  to  a  solid-like  gel 
when  exposed  to  electric  fields.  The  mechanism  for  this  transformation  is  the  formation  of 
particle  chains  across  the  field  as  displayed  in  Figure  7[l-5],  When  formed,  the  chains 
restrict  the  flow  of  the  liquid  matrix. 
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Figure  7.  ER  Phenomenon  a)no  electric  field  b)with  electric  field  applied 


Ott  et  al.  [6,7]  first  demonstrated  that  one  can  convert  chaotic  motion  to  periodic  motion. 
They  accomplished  this  by  controlling  the  system  about  one  of  the  many  unstable  periodic 
orbits  embedded  in  the  chaotic  attractor  through  small  time-dependent  perturbations  of  an 
accessible  system  parameter.  The  method,  known  as  the  OGY  algorithm,  is  different  from 
other  published  methods  that  typically  make  large  changes  in  system  behavior  and  attempt 
to  eliminate  chaos  all  together.  In  many  cases  the  latter  requires  changes  in  the  system 
that  are  too  expensive  or  impossible. 

Several  studies  have  investigated  their  techniques.  Hunt  [8]  eliminated  the  need  for  a 
Poincare  section  by  using  direct  system  measurement.  This  made  the  method  converge 
much  faster  enabling  the  manipulation  of  rapidly  changing  systems.  Shinbrot  et  al.  [9] 
devised  a  method  for  forcing  the  system  into  the  orbit  of  choice  much  faster  from  any 
arbitrary  initial  condition.  Dressier  and  Nitsche  [10]  were  able  to  improve  overall  system 
performance  by  taking  into  account  past  parameter  changes.  Ditto  et  al.  [11,12]  achieved 
the  first  experimental  control  of  chaos  by  controlling  a  buckling  magnetoelastic  beam. 
Hall  et  al.  [13]  investigated  chaotic  control  of  a  simple  buckled  beam.  They  stabilized  the 
motion  by  moving  the  system  out  of  the  chaotic  parameter  region.  Rollins  et  al.  [14] 
discussed  a  recursive  proportional  feedback  algorithm  for  highly  dissipative  systems  using 
an  adaptation  from  Dressier  and  Nitsche  [15],  Garfinkel  et  al.  [16]  successfully 
controlled,  where  previous  methods  failed,  heart  arrhythmia  in  animals  using  an  OGY 
pacemaker.  The  system  we  have  controlled  was  first  examined  by  Moon  and  Holmes 
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[17].  They  developed  a  computer  simulation  of  the  magnetoelastic  attractor  and 
compared  those  results  with  experimental  testing.  Their  results  showed  qualitative 
similarities. 

Rheological  Solutions  has  divided  the  product  line  into  discrete  and  distributed  devices. 
Discrete  devices  absorb  energy  at  a  point  source.  Dashpots,  shock  absorbers,  and  engine 
mounts  are  types  of  discrete  devices.  They  are  useful  when  the  source  of  noise  is  isolated. 
The  idea  is  to  stabilize  the  vibration  before  the  energy  reaches  its  target.  Distributed 
devices  are  used  when  we  need  to  change  the  behavior  of  an  entire  structure.  Distributed 
damping  methods  include  constrained  layer  damping  and  damping  coatings.  These 
systems  absorb  energy  from  the  structure  as  a  whole  and  are  ideal  for  attenuating  low 
amplitude  noise  at  multiple  locations  and  high  frequencies.  Their  application  is  limited  by 
their  inability  to  absorb  large  strain  amplitudes  and  are  generally  less  efficient  than  discrete 
devices.  For  this  investigation,  we  studied  the  control  of  a  distributed  Chaotic  Motion 
Stabilizer. 
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m.  RESULTS 


The  basic  recipe  for  chaotic  control  contains  three  steps.  First,  identify  periodic  orbits 
embedded  within  the  attractor.  Second,  determine  the  local  dynamics  of  those  orbits. 
Finally,  stabilize  an  orbit  around  a  fixed  point  coordinate. 

IDENTIFY  EMBEDDED  PERIODIC  ORBITS 

As  previously  discussed,  the  analysis  of  chaotic  systems  are  best  illustrated  using  a  state- 
space  representation  .  The  state  variables  appear  as  first-order  equations, 

^  =  F(X,p)  (1) 

at 

where  p  is  the  control  parameter.  A'  is  a  set  of  physical  quantities  that  completely  specifies 
the  state  of  the  system. 


*=(*o(0>*i(0,*2(0>---,*m-i(0)  (2) 

and  m  is  the  number  of  state  variables.  One  can  observe  many  unstable  periodic  orbits 
embedded  within  chaotic  attractors.  We  define  an  orbit  as  a  series  of  points  £,,£2, £,,... 
such  that, 

maps  a  Poincare  section  or  “slice”  of  the  chaotic  attractor  taken  when  one  state  variable  is 
equal  to  a  constant.  A  periodic  orbit  occurs  when  the  series  begins  to  repeat  itself.  The 
period  rth  orbit  is  defined  by, 

(4) 

DETERMINE  THE  LOCAL  DYNAMICS  OF  THE  ORBIT 

For  simplicity  we  will  examine  the  dynamics  around  a  period  one  orbit,  but  the  analysis 
also  applies  to  higher  orbits.  Let  =0  be  the  fixed  point  along  the  periodic  orbit.  By 
changing  the  control  parameter,  p ,  the  fixed  point  coordinates  will  translate  to*  some  other 
point  close  to  the  original  fixed  point,  (/?).  Assuming  that  p= 0  along  the  original  orbit 
and  that  the  changes  in  the  system  are  small, 

8  =  %f(p)  =  Zf(p)  /  P  (5) 

Using  this  approximation,  the  translation  can  then  be  linearized  about  the  original  orbit, 
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(6) 


where  M  is  an  m  x  m  matrix.  The  dynamics  around  the  point  are  found  by^substituting 
Equation  5  into  Equation  6  and  noting  that  the  M  matrix  contains  the  direction  vectors  of 
the  unstable  and  stable  manifolds, 


fc+,}  =  {R*}  +  [W-  +  -Png)  0) 

where  X u  and  Xs  are  the  unstable  and  stable  eigenvalues,  eu  and  et  are  the  unstable  and 
stable  unit  eigenvectors,  and  fu  and  f,  are  contravariant  basis  vectors  defined  by 

fueu  =f,e,  =  1  (8) 

and 

/.^=/,e«=°  (9) 


CONTROLLING  THE  SYSTEM 


Figure  8.  Chaotic  Control  of  the  nth  Iteration 


To  stabilize  the  system,  we  drive  the  orbit,  ,  towards  the  stable  manifold  on  the  desired 
fixed  orbit,  sO,  as  illustrated  in  Figure  8.  This  is  accomplished  by  manipulating  the 
accessible  control  parameter,  p,  such  that 


/.*£*  i=0 


(10) 
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The  control  parameter  required  to  accomplish  this  can  be  calculated  by  substituting 
Equation  10  into  Equation  7, 


p.  =  K(K-iY(S,f.)l(gf.)  (ID 

As  the  trajectory  approaches  this  stable  manifold,  the  system  will  collapse  towards  the 
fixed  point.  Once  on  the  manifold  we  set  p=  0. 

The  operational  range  of  such  an  algorithm  is  limited  by  the  ability  to  vary  the  control 
parameter.  Only  small  changes  in  the  control  parameter  are  desired, 


P*>P>  - p . 


(12) 


where  p-  is  the  control  threshold.  The  threshold  is  made  small  enough  that  the  system’s 
overall  dynamics  are  not  altered  by  variations  in  the  parameter  --  we  do  not  want  to  create 
new  orbits  with  different  dynamics.  Consequently,  control  is  only  enabled  when 


\tn-  Zf\<p\{1- g-  fu 


(13) 


along  the  unstable  manifold.  Outside  of  this  range  we  will  essentially  turn  the  control  off, 
recognizing  that  the  system’s  chaotic  nature  will  eventually  drive  it  back  into  the 
controllable  limits. 


HENON  A  TTRACTOR  MODEL 

The  Henon  attractor  control  provides  an  excellent  example  of  chaotic  control  methods. 
Let  x  and  y  be  coordinates  of  a  two  dimensional  Poincare  section  where, 

xn+i=A-x2n  +  By „  (14a) 

y~  i  =  (14b) 

♦ 

and  A  is  the  control  parameter  and  B  is  held  constant.  Figure  9  illustrates  the  chaotic 
nature  of  the  Henon  attractor.  We  set  a  nominal  value,  A0 ,  for  the  control  parameter  and 
vary  it  by  quantity  by  p. 


A  =  A0+p  (15) 

Our  goal  is  to  control  the  system  about  a  period  one  orbit  within  the  attractor  so  that, 

<r, (i6) 

Setting  Ao=lA  and  5=0.3,  the  period  one  orbit  results  in  a  fixed  point  located  at, 
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4r  =  =  (0884,0.884) 


(17) 


The  local  dynamics  around  this  point  are  linearized  using  Equation  6, 

fcH'fcl  *  °8) 

where  M  is  a  2  x  2  direction  matrix.  In  this  case,  the  directional  matrix  is  calculated 
analytically, 


[M].= 


A+l^n  fon+xlfyn 

-2-  x. 

0.3' 

\ldVn. 

4.=i/ 

1 

0 

(19) 


Substituting  Equation  17  into  Equation  19,  the  corresponding  unstable  and  stable 
eigenvalues  and  vectors  around  the  fixed  point  are  shown  in  Table  2  along  with  the 
eigenvalues  and  vectors  of  a  period  two  and  period  four  orbit.  We  used  the  same  analysis 
to  determine  the  local  dynamics  around  those  orbits  where, 

C:;Hm]’[H]}  (2o) 

for  period  two  orbits,  and 

M  =  [W]4[W],[M]![M],|^|  (21) 

for  period  four  orbits. 


ORBIT 

FIXED 

UNSTABLE 

STABLE 

NO. 

POINTS 

EIGEN 

VECTORS 

EIGEN 

VALUES 

EIGEN 

vector! 

EIGEN 

VALUES 

1 

(0.884,0.884) 

(-0.887,0.461) 

-1.924 

(0.163,1.043) 

0.156 

2 

(-0.666,1.366) 

(-0.928,0.373) 

-3.009 

(0.297,-1.197) 

-0.030 

4 

(0.894,-0.989) 

(0.884,0.467) 

-8.637 

(0.165,1.218) 

-0.001 

Table  2.  Unstable  and  Stable  Eigenvalues  and  Vectors  for  the  Henon  Attractor 
Governed  by  Equation  1 1,  the  size  of  the  controllable  range  for  the  period  one  orbit, 
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\Zh-4f\<  0.1084  (22) 

Applying  Equation  1 1  eveiy  time  the  state-space  trajectory  pierces  this  boundary  forces 
the  dynamical  behavior  of  the  system  to  eventually  collapse  onto  the  fixed  point.  Figure 
10  shows  the  effects  of  our  control  on  the  Henon  attractor. 


Figure  10.  Chaotic  Control  of  Henon  Attractor 
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er-magnetoelastic  attractor  model 


The  main  objective  for  Phase  I  research  is  to  simulate  a  computer  model  of  a  real  physical 
structure  incorporating  ER  materials.  We  want  to  evaluate  the  feasibility  of  introducing 
chaotic  control  algorithms  into  such  structures  for  an  experimental 
demonstration/validation  prototype.  We  chose  to  model  the  magnetoelastic  attractor. 


Moon  and  Holmes[17]  first  investigated  the  magnetoelastic  attractor.  They  modeled  the 
behavior  of  a  structure  using  a  cantilevered  Bemoulli-Euler  beam  vibrating  in  a  magnetic 
field.  Their  model  takes  into  account  the  kinetic  energy  of  the  beam,  potential  energy  of 
the  beam,  potential  energy  of  the  magnetic  field,  and  external  work.  For  our  project,  we 
have  modified  their  equations  to  account  for  ER  material  adaptability  and  an  external 
electromagnetic  forcing  function. 


ER  material  adaptability  manifests  itself  in  the  flexural  rigidity  of  the  structural  system. 
We  modeled  this  parameter  based  on  a  composite  structure  with  ER  material  sandwiched 
between  two  elastic  constraining  layers.  Don  et  al.[18-20]  investigated  the  theoretical 
dynamic  behavior  of  such  structures  using  the  Ross,  Kerwin,  and  Ungar[21]  constrained 
layer  damping  models.  The  model  accounts  for  ER  material  properties  by  assuming 
shearing  motion  within  the  sandwiched  layer  and  incorporating  these  properties  into  an 
effective  flexural  rigidity.  Equation  23  is  the  effective  flexural  rigidity, 


where 


El  =  £,/,  +  E2I2  +  £3/3  +  EXHXH?0  +  E2H2H]0  +  E3H3Hj0 


e2i2h3X-d 

H2  1  +  g 


F  TT 

I-J2*12  Tj  1  tr  JJ  IT 
**20~*~  £'3*1 3*1' 


30 


Hn-D 


l  +  g 


(23) 


D  E2H2{H2X  -  H3X  /  2)  +  g(E2H2H2X  +  E3H3H3X) 
EXHX  +  E2H2  /  2  +  g(ExHx  +  E2H2  +  E3H3) 


(24) 


and 


e3h3h2£ 


(25) 


£/,  is  the  flexural  rigidity  of  the  rth  layer,  X  is  the  wave  number,  and  H\0  is  the  distance 

from  the  mid  plane  of  the  /th  layer  to  a  new  neutral  plane  created  by  the  addition  layers 
two  and  three.  The  old  neutral  plane,  defined  as  the  neutral  plane  without  any  additional 
layers,  is  the  mid  plane  of  the  first  layer.  Therefore, 

Hi0  -  Hn-D  (26) 
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where  Hu  is  the  distance  between  the  mid-plane  of  the  /th  layer  and  the  mid-plane  of  layer 
one,  and  D  is  the  displacement  of  the  old  neutral  plane.  Equation  25  incorporates  shear 
behavior  of  the  sandwiched  layer  where  the  value  is  a  function  of  both  the  ER  material  and 
a  sealant  used  to  hold  the  fluid, 

G2  =  ^f*--G„^+-£sL-Gm  (27) 

™1otaI  ™ total 

where  A  represents  the  surface  area  of  the  respective  material.  The  shear  modulus  of  the 

ER  material  is  found  experimentally  to  be  a  function  of  the  electric  field  squared, 

* 

Ger  =  K-V2  (28) 


where  K  is  an  experimentally  determined  material  constant  and  V  is  the  applied  electric 
field. 


■ 

I 

mMmmm 

DESCRIPTION 

VALUE 

K 

ER  Material  Property  Constant 

15.34  kPa/kV/mm 

Ea 

Elastic  Modulus  of  ER  material 

OPa 

E 

Elastic  modulus  of  layers  1  and  3 

2.34  GPa 

Gj taUmt 

Shear  modulus  of  the  sealant  material 

0.5  MPa 

H, 

Height  of  the  1st  layer 

7.94  x  10-4  m 

h2 

Height  of  the  2nd  layer 

7.94  x  10-4  m 

Hs 

Height  of  the  3rd  layer 

7.94  x  10^  m 

k 

Wavenumber 

10.29 

w 

Width  of  the  structure 

0.025  m 

L 

Length  of  the  structure 

0.305  m 

V 

Structural  damping  per  unit  length 

.038  kg/s/m 

P 

Magnetic  field  constant 

0.5  kg/m’/s1 

a 

Magnetic  field  constant 

18.32  kg/m/s1 

Ao 

External  Force  Amplitude 

0.28TN 

n 

Forcing  frequency 

0.89  rad/s 

m 

Mass  per  unit  length 

0.115  kg/m 

Table  3.  Physical  Constants  for  the  ER-Magnetoelastic  Model 


We  substituted  this  effective  flexural  rigidity  into  a  modified  Moon  and  Holmes 
magnetoelastic  system.  The  Moon  and  Holmes  model  assumed  a  moving  reference  frame 
as  the  external  forcing  function.  For  this  investigation,  a  discrete  forcing  function  actuates 
the  structure  at  a  point  source.  The  resulting  equation  of  motion  for  the  system  is, 

m(x)w  +  tjw-ocw  +  flw3  +  EIw  =  A0  sin(Q/)  (29) 

where  m(x)  is  the  mass  per  unit  length,  w  is  the  horizontal  position,  ij  is  the  overall 
structural  damping,  A0  is  the  magnitude  of  the  forcing  function,  and  a  and  /?  are  magnetic 
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field  constants.  This  model  assumes  only  first  mode  vibrational  behavior  and  magnetic 
field  affects  only  at  the  beam’s  tip.  The  state  equations  representing  the  dynamical 
equation  of  motion  in  Equation  29  are, 

ti'o  =  (-rpv0  +  awx  -  pw\  -  EIwx  +  w2 )  /  rri(x) 

wx  =  w0  (30a, b,c) 

w2  =  ^40Qcos(ifl/) 

where  wa  is  the  velocity  of  the  beam,  w,  is  the  position  of  the  beam,  and  w,  the  external 
forcing  function.  The  physical  modeling  constants  are  defined  in  Table  3. 

A  three  dimensional  attractor  shown  in  Figure  1 1  is  found  by  iterating  Equation  31, 

W(t  +  A t)  =  f(W(t) )  •  At  +  Wit)  (31) 


where  W  is  the  set  of  state  variables.  Figure  12  shows  a  surface  section  or  Poincare 
section  of  this  three  dimensional  attractor  taken  at 


2701 


a 


(32) 


A  period  three  orbit  is  observed  from  this  mapping  such  that, 

£,=Cb  (33) 

We  selected  to  control  the  system  using  the  fixed  point  located  at, 

&  =  (Wf)  =  (-0.404,-0.073)  (34) 


The  corresponding  eigenvalues  and  vectors  are  found  by  evaluating  the  behavior  of  several 
points  around  the  fixed  point.  Table  4  describes  the  local  dynamics. 


1  UNSTABLE 

EIGENVECTORS 

EIGENVALUES 

EIGENVECTORS 

EIGENVALUES 

(0.284,0.959) 

23.00 

(-0.595,0.804) 

0.13 

Table  4.  Unstable  and  Stable  Eigenvalues  and  Vectors  for  the  Magnetoelastic  Attractor 


Figure  13  and  14  show  the  effects  of  our  control  algorithm  on  the  state-space  and  time 
domains  respectively. 


18 


Figure  12.  Computer  Simulated  Poincare  Section  of  the  Magnetoelastic  Attractor 
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EXPERIMENTAL  MODEL 


Our  feasibility  analysis  includes  a  comparison  between  this  computer  simulation  and  an 
experimental  model.  We  fabricated  an  electrorheological  material  adaptive  structure 
composed  of  Lord  Corporation  Material,  ER  in,  sandwiched  between  two  polycarbonate 
sheets.  Neoprene  rubber  adhered  to  the  edges  of  the  polycarbonate  sealed  the  structure 
and  acted  as  a  spacer  between  the  constraining  layers.  The  adhesive  was  a  3M  847 
Rubber  and  Gasket  adhesive.  A  graphite  coating,  Grapho  206,  applied  to  the  surfaces  of 
the  polycarbonate  sheets  served  as  the  electrodes.  Figure  15  shows  a  schematic 
illustration  of  the  structure  and  its  dimensions. 

A  detailed  illustration  of  the  experimental  testbed  is  shown  in  Figure  16.  A  Bently  Nevada 
3040  HTB  electromagnetic  probe  actuated  the  structure  while  an  AROMAT  30mm  Laser 
Proximity  Sensor  measured  the  vibrational  amplitude.  Two  rectangular  magnets  attached 
on  a  platform  at  the  end  of  the  beam  produced  the  magnetic  field.  All  relevant  data  was 
recorded  using  a  PC  based  acquisition  system  supplied  by  National  Instruments. 

The  Figures  17  and  18  illustrate  the  chaotic  behavior  and  Poincare  Section  Respectively, 
while  Figure  19  shows  the  orbits  we  desire  to  stabilize. 


Figure  15.  Experimental  ER  Material  Composite  Structure 
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FUNCTION  GENERATOR 


ELECTROMAGNETIC 

ACTUATOR 


Figure  16.  Schematic  of  the  Experimental  Testbed 
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IV.  SUMMARY  OF  IMPORTANT  RESULTS 


•  Periods  one,  two,  and  four  orbits  within  a  Henon  attractor  were  effectively  stabilized. 

•  The  chaotic  algorithm  applied  to  the  ER-magnetoelastic  computer  model  effectively 
stabilized  the  system  about  a  period  three  orbit. 

•  The  computer  model  of  the  ER-magnetoelastic  system  agreed  qualitatively  with  the 
experimental  simulation. 

•  We  observed  period  three  and  period  five  orbits  within  the  experimental  ER- 
magnetoelastic  attractor 


V.  DISCUSSION 

The  most  important  result  from  our  research  is  that  chaotic  controls  can  be  implemented 
onto  a  computer  simulated  model  and  that  the  simulation  qualitatively  resembles  a  real 
physical  structure.  Chaos  is  a  “visual”  science  where  qualitative  behavior  is  more 
important  than  quantitative  equivalence.  Exact  “scale”  between  theory  and  experiment  is 
trivial.  These  simulations  are  used  as  tools  to  quantify  qualitative  behavior  where  we  ask  - 
-  What  degree  of  chaos  and  types  of  attractors  can  we  control? 

Though  there  are  visual  differences  in  the  attractor  models,  the  forms  of  the  attractor  are 
qualitatively  similar.  A  comparison  of  the  theoretical  and  experimental  attractors  shown  in 
Figures  11  and  17  illustrate  this  similarity.  The  noticeable  differences  most  likely  result 
from  the  experimental  structure’s  deviations  from  pure  sinusoidal  mode  shapes  and  the 
reliability  of  the  rheological  property  data.  The  theoretical  model  uses  an  approximation 
of  an  effective  flexural  rigidity  of  a  3  layer  composite  structure.  The  flexural  rigidity  is 
based  on  the  Ross,  Kerwin,  and  Ungar  (RKU)  [21]  damping  model  that  only  assumes  pure 
sinusoidal  mode  shapes.  A  cantilevered  beam  does  not  have  pure  sinusoidal  mode  shapes 
especially  at  lower  frequencies.  In  this  investigation,  we  are  actuating  a  cantilevered 
structure  at  frequencies  where  mode  1  behavior  dominates.  We  believe  some  of  the 
differences  in  theoretical  and  experimental  results  arise  from  this  RKU  assumption.  In 
addition,  Yalcintas  et  al.  [20]  have  questioned  the  reliability  of  the  rheological  behavior 
data  obtained  on  these  ER  materials.  More  rheology  testing  is  clearly  needed. . 

Since  we  are  able  to  control  the  computer  model  of  the  ER-magnetoelastic  system  and  the 
experimental  model  qualitatively  agrees  with  the  computer  simulation,  we  believe  that 
those  same  chaotic  algorithms  can  control  the  physical  experiment.  Figures  14  and  19 
show  the  theoretical  and  experimental  orbits  we  desire  to  control. 


VL  CONCLUSION 

We  believe  we  can  stabilize  the  chaotic  attractor  of  the  ER  magnetoelastic  system. 
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VH.  FURTHER  RESEARCH  AND  DEVELOPMENT 

There  are  several  issues  that  need  to  be  addressed  before  we  can  bring  products  like  this 
into  the  marketplace.  The  following  section  defines  our  objectives  and  the  questions  we 
need  to  answer  before  a  commercial  product  could  become  available. 

1.  Simulate  Chaotic  Vibration  Control 

•  What  types  of  chaotic  attractors  can  we  control? 

•  What  are  the  best  algorithms  to  use? 

•  What  are  the  limitations  of  chaotic  control? 

2.  Fabricate  Prototype  Chaotic  Motion  Stabilizers 

•  What  degrees  of  chaos  can  we  control? 

•  What  are  the  physical  limitations? 

•  How  robust  is  the  controller? 

•  How  will  these  systems  be  manufactured? 

3.  Integrate  Stabilizer  with  Controls  for  Real  World  Applications 

•  Are  chaotic  vibrations  effective  in  real  world  environments? 

4.  Investigate  the  Rheological  Behavior  of  ER  Materials 

•  What  are  the  best  fluids  to  use  for  each  application? 

•  How  reliable  is  the  long  term  fluid  behavior? 

•  What  are  the  effects  of  temperature  on  fluid  behavior  and  reliability? 

5.  Investigate  other  Methods  to  Change  the  System  Parameters 

•  What  other  controllable  materials  could  we  use  to  implement  a  Chaotic  Motion 
Stabilizer? 


27 


Vni.  REFERENCES 


1.  Block,  H.  and  J.P.  Kelly,  Electro-rheology.  J.  Phys.  D:  Appl.  Phys.,  1988. 1:  p.  1661-1677. 

2.  Gast,  A.P.  and  C.F.  Zukoski,  Electrorheological  Fluids  as  Colloidal  Suspensions.  Advances  in 
Colloid  and  Interface  Science,  1989.  30:  p.  153-202. 

3.  Jordan,  T.C.  and  M.T.  Shaw,  Electrorheology.  IEEE  Transactions  on  Electrical  Insulation,  1989. 
24(5):  p.  849-878. 

4.  Weiss,  K.D.,  J.P.  Coulter,  and  J.D.  Carlson.  Electrorheological  Materials  and  Their  Usage  in 
Intelligent  Material  Systems  and  Structures,  Parti:  Mechanisms,  Formulations  and  Properties,  in  Recent 
Advances  in  Adaptive  and  Sensory  Materials  and  Their  Applications.  1992.  Blacksburg,  Virginia: 
Technomic  Publishing  Co.,  Inc. 

5.  Coulter,  J.P.,  K.D.  Weiss,  and  J.D.  Carlson.  Electrorheological  Materials  and  Their  Usage  in 
Intelligent  Material  Systems  and  Structures,  Part  II:  Applications,  in  Recent  Advances  in  Adaptive  and 
Sensory  Materials  and  Their  Applications.  1992.  Blacksburg,  Virginia:  Technomic  Publishing  Co.,  Inc. 

6.  Ott,  E.,  C.  Gregobi,  and  J.A.  Yorke,  Controlling  Chaos.  Physical  Review  Letters.  Physical  Review 
Letters,  1993.  64(11):  p.  1196. 

7.  Ott,  Grebogi,  C.,  Romeiras,  F.,  and  Yorke,  J.A.,  Critical  Exponents  for  Crisis-Induced  Intermittancy, 
Physical  Review  letters,  1987,  36(1  l):p.  5365. 

8.  Hunt,  E.R.,  High-Period  Orbits  Stabilized  in  the  Diode  Resonator,  1st  Experimental  Chaos  Conference, 
Arlington,  VA,  1991. 

9.  Shinbrot,  T.,  Grebogi,  C.,  Ott,  E.,  Yorke,  J.A.,  Using  Small  Perturbations  to  Control  Chaos,  Nature, 
1993,  363(6428):  p.411-17. 

10.  Dressier,  U.  and  Nitsche,  G.,  Controlling  Chaos  Using  Time  Delay  Coordinates,  Physical  Review 
Letters,  1991,  86(1):  p.l. 

11.  Ditto,  W.L.,  S.N.  Rauseo,  and  M.L.  Spano,  Mastering  Chaos.  Scientific  American.  Scientific 
American,  1990.  (August):  p.  78. 

12.  Ditto,  W.L.,  S.N.  Rauseo,  and  M.L.  Spano,  Experimental  Control  of  Chaos.  Physical  Review  Letters. 
Physical  Review  Letters,  1990.  65(26):  p.  3211. 

13.  Hall,  E.K.,  and  Hanagud,  S.V.,  Control  of  Nonlinear  Structural  Dynamic  Systems:  Chaotic 
Vibrations,  Journal  of  Guidance,  Control,  and  Dynamics,  16(3):  P.  470-6. 

14.  Rollins,  R.W.,  Parmananda,  P.,  and  Sherard,  P.,  Controlling  Chaos  in  Highly  Dissipative  Systems:  A 
Simple  Recursive  Algorithm,  Physical  Review  E,  47(2):  p.  R780-3. 

16.  Garfinkel,  A.,  Spano,  M.L.,  Ditto,  W.L.,  and  Weiss,  J.N.,  Controlling  Cardiac  Chaos,  Science,  1992, 
257:  p  1230-35. 

17.  Moon,  F.C.,  and  Holmes,  P.J.,  A  Magnetoelastic  Strange  Attractor,  Journal  of  Sound  and  Vibration, 
1979,  65(2):  p.  275-96. 

18.  Don,  D.L.,  An  Investigation  of  ER  Material  Adaptive  Structures.  1993,  Lehigh: 

19.  Don,  D.L.  and  J.P.  Coulter,  An  Analytical  and  Experimental  Investigation  of  ER  Material  Based 
Adaptive  Beams.  Journal  of  Intelligent  Materials  and  Systems.  Journal  of  Intelligent  Materials  and 
Systems,  Submitted  1993. 

20.  Yalcintas,  M.,  J.P.  Coulter,  and  D.L.  Don,  Structural  Modeling  and  Optimal  Control  of  ER  Material 
Based  Plates.  ASME  Journal  of  Vibration  and  Acoustics.  ASME  Journal  of  Vibration  and  Acoustics, 
Submitted  1993. 

21.  Ross,  D.,  E.E.  Ungar,  and  E.M.  Kerwin,  Damping  of  Plate  Flexural  Vibrations 

by  means  of  Viscoelastic  Laminate.  Structural  Damping.  Structural  Damping,  1959.  ASME:  p.  49-88. 

22.  Crutchfield,  J.P.,  Farmer,  J.D.,  Packard,  N.H.,  Shaw,  R.S.,  Chaos,  Scientific  American,  1986,  255: 
p46-57. 


28 


